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Collective excitations of a dipolar Bose-Einstein condensate 
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We have measured the effect of dipole-dipole interactions on the frequency of a collective mode of 
a Bose-Einstein condensate. At relatively large numbers of atoms, the experimental measurements 
are in good agreement with zero temperature theoretical predictions based on the Thomas Fermi 
approach. Experimental results obtained for the dipolar shift of a collective mode show a larger 
dependency to both the trap geometry and the atom number than the ones obtained when measuring 
the modification of the condensate aspect ratio due to dipolar forces. These findings are in good 
agreement with simulations based on a gaussian ansatz. 

PACS numbers: 03.75.-b, 03.75.Kk, 67.85.De 
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Interactions strongly affect the properties of quantum 
degenerate gases [l| • While in most of the experiments to 
date short range interactions dominate, the recent pro- 
duction of Bose-Einstein condensates (BEC) of highly 
magnetic Chromium atoms [2j, |3j, and the tremendous 
rogresses in the manipulation of hetcronuclear molecules 
have brought considerable interest towards particles 
interacting through long-range, anisotropic, dipole-dipole 
interactions (DDI) [5]. 

The analysis of collective excitations is an excellent 
tool to analyze the effects of interactions in many-body 
systems [l|, |6j. In trapped BECs where short range 
isotropic interactions dominate, the excitations at low en- 
ergies, being analogous to phonons in homogeneous sys- 
tems, have a collective character. These excitations are 
well understood in the framework of the Gross-Pitaevskii 
equation (non linear Schrodinger equation). DDI add a 
non-local character to this non-linear framework. In 3D, 
it has been predicted that relatively small non-local in- 
teractions from dipoles should lead to modifications of 
collective modes frequencies. This effect of DDI, which 
depends on the orientation of the polarization axis of the 
dipoles with respect to the trap geometry [7|, |8[ , has to 
date not been observed. 

For larger intensities of DDI in 3D, or in reduced di- 
mensions, the effect of DDI on the properties of BECs is 
dramatic. In 3D, both the ground state and the collec- 
tive excitations are strongly modified 9] , eventually lead- 



ing to collapse IIJ lllj. In 2D, the excitation spectrum 
presents a roton structure reminiscent of the physics of 
liquid helium [12j . This roton minimum vanishes in ID, 
leading to the breakdown of the Landau criterion for su- 
perfluidity Il3| .|l4l| (similar to what is found with contact 
interactions [15|). Understanding the collective modes of 
trapped dipolar quantum gases has a peculiar interest be- 
cause both collective modes and long range interactions 
are key components in the quantum computing toolbox 

The recent productions of BECs with chromium atoms 
[2j, y] , which have a relatively large magnetic dipole mo- 



ment (6 Bohr magnetons), makes it possible to experi- 
mentally investigate the effects of DDI on the properties 
of quantum degenerate gases. Up to now, the only mea- 
sured effect of DDI on the hydrodynamic properties of 
BECs has been the modification of the Thomas-Fermi 
radii of the BEC [H, and its implosion if DDI be- 
come comparable to short range interactions [lfj . In this 
Letter, we report the first measurement of the modifica- 
tion of collective excitations frequencies due to DDI. 

In order to experimentally show that a shift of a low 
energy collective excitation frequency is induced by DDI, 
we have measured the frequency of a collective mode of 
our spin polarized Cr BEC for two orthogonal orienta- 
tions of the magnetic dipoles. The measured shift is very 
sensitive to trap geometry, a consequence of an interplay 
between the anisotropies of the trap and of DDI. De- 
spite the fact that the number of atoms in the BEC is 
relatively small, the observed shift is in good agreement 
with zero temperature theoretical predictions based on 
the Thomas Fermi (TF) approximation 8] . By operating 
at even lower number of atoms, we measure a shift which 
decreases, and the deviation from the TF prediction is 
faster than what is observed for the aspect ratio of the 
expanded BEC, in good agreement with our numerical 
simulations. 

Our all-optical method to produce Chromium BECs is 
described in 3] . Its main specificity is the direct loading 
of an optical trap from a Magneto Optical Trap. Atoms 
are first loaded in a horizontal (x axis) infra-red (IR) one- 
beam optical dipole trap produced by a 1075 nm, 50 W 
fiber laser. Then some of the IR light is transferred to 
a vertical beam (z axis) by rotation of the angle of a 
A/2 wave plate, to produce a strongly confining trap in all 
three directions of space, where evaporation is performed. 
The horizontal and vertical beam waists are respectively 
equal to 40 and 50 fim, and 4> — when all the laser power 
(30 Watt) is in the horizontal beam. In contrast to |3j, 
the horizontal trapping beam is not retro-reflected, which 
increases the stability of the trap: the r.m.s. fluctuations 
of the BEC TF radii are reduced from typically 10 per- 
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FIG. 1: Influence of DDI on the intermediate collective mode frequency. Free oscillations of the mode after an excitation in a trap set by <fi — 27° 
is monitored by measuring the aspect ratio of the BEC, given by the two experimental TF radii along the y and z axis. The magnetic field is cither 
vertical (black diamonds), or parallel to the horizontal imaging beam (red dots). Wc plot as well the best fits from damped sinusoidal forms. The 
residuals to the fit arc shown below. Inset: schematic of the experiment. 



cent to below 3 percent. Condensates can be produced 
in different trap geometries characterized by oscillation 
frequencies w XiJ/iZ , by adiabatically changing <j) after the 
quantum regime is reached. 

Once the BEC is obtained in a trap for a given value of 
(/>, we excite shape oscillations of the BEC by modulating 
the IR laser power with a 20 % amplitude during about 
15 ms, at a frequency close to the second surface collec- 
tive mode resonance (here referred to as the intermediate 
collective mode). This quadrupole-like mode corresponds 
to an oscillation out of phase of the different TF radii 
(see fig. [T]). We then let the cloud oscillate freely in 
the crossed optical trap during a variable time. We fi- 
nally switch off the trap, let the BEC expand for 5 ms, 
and take destructive absorption images thus measuring 
the values of the cloud TF radii. We only report results 
for the intermediate collective mode: the lower collective 
mode is predicted to have a larger shift due to DDI, but 
we found that for our trap parameters this mode is very 
hard to excite selectively, due to a quasi-degeneracy be- 
tween its second harmonic and the higher surface mode 
(monopole like) frequency [21]. 

A key feature of dipole interactions is its anisotropic 
nature. In an anisotropic trap, DDI therefore depends 
on the orientation of the spins (set by a magnetic field) 
relative to the axis of the trap. We use this property 
to perform a differential measurement of the shift of the 
intermediate collective mode due to DDI: we measure 
the oscillations of the condensate aspect ratio for two 
perpendicular orientations of the magnetic field 6. For 
6 = 0, the B field is horizontal and quasi aligned with 
both the imaging beam and the IR trapping beam, while 
for 6 = ir/2 it is vertical. We then fit the curves by 
exponentially damped sinusoidal functions, and deduce 



the corresponding oscillation frequencies ujq(8). Both 
the experimental results and the fits are shown in fig [T] 
We also show the residuals of the fit. The noise on the 
experimental data corresponds to a typically 3 % r.m.s. 
noise on the measured TF radii. This noise is neither 
related to the fluctuations in atomic number, nor to the 
way collective excitations are produced. We stress that 
fig. [1] shows that the residuals to the fit do not increase 
with oscillation time t, a signature of very good short 
term (20 ms: the typical oscillation duration) and long 
term (20 s: the cycling time of the experiment) stability. 

As shown in fig. [TJ damping of the collective excita- 
tions is rather strong, likely due to the large anharmonic- 
ity characteristic of an optical dipole trap. The damping 
rate depends neither on the magnetic field orientation 8 
nor on the trap anisotropy set by <f>. 

To deduce the shift due to DDI on the collective mode 
frequency from the experimental data, it is important 
to measure and understand all other systematic shifts 
associated with varying 8. For this, it is in theory neces- 
sary to measure the systematic shift of all three orthog- 
onal vibrational frequencies of the trap. In practice, the 
knowledge of the shift of the vertical axis is the most im- 
portant, and we have verified that the effect of the shifts 
of the frequency of the other two modes are negligible. 
We therefore excite the dipole mode along the vertical di- 
rection, and measure the frequency of the center of mass 
oscillation of the condensate U) z (d). We deduce the rel- 
ative shift 5 D = 2(w*(0) -u z (n/2))/(u z (0) + oj z (tt/2)), 
which is plotted in fig. [2] 

Figure [2] shows that systematic effects on the verti- 
cal dipole frequency strongly depend on the trap geom- 
etry. We have found two independent sources for this 
systematic shift. The first is the presence of magnetic 
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FIG. 2: Variation with the magnetic field orientation of the 
vertical dipole mode frequency, as a function of the trap ge- 
ometry. The experimental relative shifts 8n are plotted versus 
the angle <j> (see text). A fit assuming a (constant) tensorial 
shift and a (trap geometry dependent) magnetic gradient ori- 
gin for the shifts is also shown. 



field gradients depending on the applied magnetic fields. 
In a parabolic trap, a potential gradient merely shifts 
the position of the center of the trap; in a gaussian 
trap, in addition to this spatial shift, the oscillation fre- 
quency is also modified. The relative frequency shift is 
Slu/uj — —2>g/ra 2 w 2 Lo A , where g is the acceleration due 
to the potential gradient, m the mass of the atoms, w 
the waist of the gaussian potential, and u> the trap fre- 
quency at the bottom of the trap. As lo depends on the 
trap geometry through the angle <fi, this first systematic 
shift also depends on <\>. The second source of systematic 
shift, A, is related to the tensorial light shift in Cr: even 
though the IR laser is very far detuned compared to the 
fine structure of the first electronically excited state, the 
AC stark-shift of Cr atoms slightly depends on both the 
laser polarization orientation, and on the internal Zee- 
man state m |22j. This second systematic shift is inde- 
pendent of <p. As shown in fig. [2] we fit our experimental 
data to a functional form A— a/u; 4 , and deduce the tenso- 
rial light shift of chromium for our experimental parame- 
ters A = 3 ± 1 percent, to be compared to our calculation 
A = 1.2 percent [22J. To our knowledge, this represents 
the first measurement of the tensorial light shift of Cr. 

Once we have measured the experimental shift of 
the collective mode frequency S exp = 2(loq(0) — 
wq(tt /2)) / (ujq(0) + ujq(tt/2)), we estimate the shift of 
the collective mode due to DDI, Sq, by subtracting from 



J exp 



the systematic shift of the intermediate collective 
mode frequency due to Sd (which we estimate using the 
theory of [23J ) . As shown in fig [3] we find for various 
trap geometries, a good agreement between the experi- 
mentally measured Sq, and a numerical model general- 
izing the theoretical results of [8| to non axisymmetric 
parabolic traps in the TF regime. 
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FIG. 3: Effect of the DDI on the experimental relative shifts 
induced by a 90 degree B field rotation, as a function of the 
trap geometry: 5q (red squares) is the shift of the interme- 
diate collective mode frequency, while S a (black dots) is the 
shift of the expanded BEC aspect ratio. Results of numerical 
simulations are also shown. 



We also plot in fig. [3] the measured modification of 
the aspect ratio due to DDI for the BEC measured 5 
ms after its release from the trap (similar to results re- 
ported in [20|). For this, we use the BEC shape oscilla- 
tion experimental data (see figQ]), and deduce the BEC 
equilibrium aspect ratio, <j(6), from the mean value of 
the TF radii. The corresponding relative variation is 
S a = 2(cr(0) - cr(7r/2))/(cr(0) + cr(7r/2)). We see from 
fig. [3] that 5 a is almost constant for various trap ge- 
ometries, which is in reasonable agreement with theory. 
In contrast, the shift of the collective mode Sq strongly 
depends on geometry. 

The large sensitivity of Sq to trap geometry comes from 
the fact that the sign of the shift of the collective mode 
is approximately set by the sign of the mean-field due 
to DDI at the center of the BEC 8] . As the sign of the 
meanfield due to DDI changes (in cylindrical traps) when 
the trap is modified from oblate (disk-like) to elongated 
(cigar- like), the shift of the collective mode changes sign 
too. For our non-axisymmetric trap, we measure Sq in 
a domain where the sign of dipole-dipole meanfield flips 
(A = uj x /uj z close to 1, see fig. |3])), hence the very large 
sensitivity illustrated by fig. [3] In contrast, as DDI al- 
ways stretch the BEC along the direction set by B, S a 
keeps the same sign whatever the trap geometry is, hence 
its lower relative variation with </>. 

Although our maximal number of atoms is rather 
small, our results for Sq coincide with theoretical pre- 
dictions for the shift of collective excitations due to DDI 
based on the TF approximation. In our experiment the 
mean-field due to DDI is only approximately equal to the 
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FIG. 4: Variation with the BEC atom number of the relative 
shifts induced by a 90° B field rotation, for two different trap 
geometries (blue circles: <f> — 27°; red squares: (j) ~ 26°). 
(a) variation of Sq and comparison with simulation (for cj> — 
26°). (b) variation of 8 a , with the corresponding results of our 
simulations. The experimental intermediate collective mode 
frequency is also plotted (black diamonds), the black line is a 
guide for the eye. 



quantum kinetic energy Ti /mR^ F « 25 Hz: the fact that 
our experimental results follow predictions based on the 
TF approximation is therefore not obvious. To deepen 
our understanding, we therefore repeated our measure- 
ments for BECs with even lower numbers of atoms, ob- 
tained by loading less atoms in the optical dipole trap 
before evaporation. As shown in fig. |U(a), we observe a 
rapid decrease of Sq as the number of atoms is lowered, 
marking a clear departure from the TF predictions. On 
the contrary, fig. |4] (b) shows that the shift of the aspect 
ratio due to DDI, S a , is quite insensitive to the num- 
ber of atoms. We have checked that for the variation in 
number of atoms that we have explored, the collective 
frequencies themselves show little shift compared to the 
TF predictions without DDI (see fig. 0] (b)). 

Our results show that measuring collective excitations 
is a much more sensitive probe of DDI than measuring 
the stretching of the BEC along the axis of the dipoles. 
We have performed numerical simulations based on a 
gaussian ansatz which takes into account the quantum 
kinetic energy [2J|: they confirm that Sq is more sen- 
sitive to a reduction of the number of particles than Jo- 
is. Typically, our numerical results show that it requires 
about three times more atoms to reach the TF predic- 
tions when measuring shift of collective excitations com- 
pared to when measuring stretching of the condensates. 
As can be shown in fig. [4] our numerical results are in 
relatively good agreement with experimental data. The 
slight disagreement between theory and experiment for 
the collective excitation frequency shift may indicate the 
limits of a simple gaussian ansatz. 

In conclusion, we have characterized the effect of DDI 
on a collective mode of a Cr BEC. For large enough num- 



ber of particles in the BEC, our results are well explained 
by TF predictions. In particular, we find a very large sen- 
sitivity of the collective mode shift as a function of the 
anisotropy of the trap, a consequence of the anisotropic 
character of DDI. We also find that our results signifi- 
cantly depart from TF predictions for lower numbers of 
atoms, even when the striction of the BEC due to DDI 
is still very well accounted for by a TF theory. This 
surprising feature is another example of the usefulness 
of collective modes to characterize quantum degenerate 
gases. Finally, we have measured for the first time the 
tensorial light shift of Cr atoms. 
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